Different types of integrability and their relation to decoherence in central spin models 
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We investigate the relation between integrability and decoherence in central spin models with more 
than one central spin. We show that there is a transition between integrability ensured by Bethe 
ansatz and integrability ensured by complete sets of commuting operators. This has a significant 
impact on the decoherence properties of the system, suggesting that it is not necessarily integrability 
or non-integrability which is related to decoherence, but rather its type or a change from integrability 
to non-integrability. 
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The Liouville-Arnol'd theorem states that if a system 
with n degrees of freedom has n involutive integrals of 
motion, which are functionally independent, its Hamilto- 
nian equations of motion are solvable via quadratures [l] . 
Such a system is called integrable. Despite for huge effort, 
so far it has not been achieved to adapt the concept of 
integrability to the quantum mechanical framework sat- 
isfactorily. At the present time there are two commonly 
accepted definitions: A quantum mechanical system is 
called integrable (i) if there is a Bethe ansatz [2 or (ii) 
if the system has a complete set of commuting operators 
(CSCO) [3 sharing "suitable" properties (to be further 
explained below) . Note that the notion of integrability in 
classical mechanics does not require the solvability of the 
quadratures. In this sense both of the aforementioned 
approaches are in direct analogy with classical mechan- 
ics. 

In investigations mainly focused on the first type of 
integrability, evidence has been found that it is related 
to transport properties to quantum phase transitions 
[5] , and to decoherence ^ [7] . Here systems of the form 

H — He + Hc^h + • • • further terms (1) 

have been considered, where He denotes a central system 
and Hc^i) a coupling term between the central system 
and a bath. Mainly two roads have been followed. On 
the one hand, the influence of chaotic or regular baths on 
the decoherence of the central system has been investi- 
gated |6 . On the other hand, the decoherence properties 
of the central systems of models which are integrable or 
non-integrable have been studied [7]. The usual proce- 
dure within such considerations is to evaluate numerically 
the level statistics of the respective system and to relate 
a possible change in the statistics to a change of other 
properties of the system happening at the same point. 

Motivated by their important role in the context of 
solid state quantum information processing [8] , we inves- 
tigate in the present letter integrability and its relation 
to decoherence in central spin models. Here we define 
a quantum system to be integrable if it is possible to 
compute all eigenstates and eigenvalues of the respec- 
tive Hamiltonian using operations with less complexity 



than the direct diagonalization of the Hamiltonian ma- 
trix [9 . Here we refer to the computional complexity. 
The exact diagonalization of a Hamiltonian matrix for 
example grows exponentially with the system size. This 
very strict notion of integrability contains (i) and (ii) as 
possible sources of integrability. 

First we study the integrable structure of central spin 
models. In particular we show that there is a transition 
between integrability ensured by Bethe ansatz and inte- 
grability ensured by CSCO. Differently from the previous 
investigations described above, we then open a new route 
by applying a strong magnetic field to the central spin 
system, and analyze its reaction with respect to deco- 
herence. In the non-integrable case as well as in the case 
of integrability ensured by Bethe ansatz, the strong mag- 
netic field leads, as generally expected, to highly coherent 
central spin dynamics, whereas in the remaining case de- 
coherence still takes place. In contrast to previous work 
we relate the latter observation explicitly to the type of 
integrability and interpret the result from two different 
points of view. 

The Hamiltonian of a central spin model is given by 

N 

i=l j=i-\-l k — 1 

+ ^^^Jij^i'^j^ (2) 

i<j 

where in the following we consider Jij = J and A^c > 1- 
For later convenience we define A = ^k=i ^j=i ^k- 
In the second identity we rewrote the original Hamilto- 
nian into terms of sums and differences between the dif- 
ferent central spins. The first term is nothing else than a 
Gaudin model [10 with a central spin replaced by a sum 
over a set of spins, whereas the second term acts as a per- 
turbation, vanishing whenever A\ = A^. =: A^. Hence it 
has to be expected that this case is integrable, whereas 



2 



the model generally should be non-integrable. This pre- 
diction has been verified explicitly in [TT by a detailed 
investigation of the spectral statistics of the model. We 
will come back to the integrable case of two central spins 
with A} = =: Ai below. 

Let us first, however, investigate in more detail general 
features of the above system, fulfilling = A^. =: A^. 
The central spins can couple to different values of the 

total central spin squared S'^ = (^^^i Si^ . Fixing 

the associated quantum number S and defining |0) = 
\S)\Ii, . . . , /at), we arrive at a usual Gaudin model with 
eigenstates 



i=l \ 3 = 1 J 



(3) 



and eigenvalues 



^=l 3=\ 

The parameters cj^ are determined by the Bethe ansatz 
equations: 
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Here Nd is the number of spin fiips compared to |0) [12]. 
Note that these equations are valid for any spin length S 
and hence any number of central spins Nc. Considering 
the Bethe ansatz equations instead of the direct diago- 
nalization of the Hamiltonian matrix reduces a problem 
of exponential complexity to one of polynomial complex- 
ity. Hence the Hamiltonian ^ with A\ = A^^ =\ A^ is 
integrable, provided the Bethe ansatz equations yield the 
correct number of solutions {cji, . . . ,co'^}. This however 
strongly depends on the inhomogeneity of the couplings 
Ak. Indeed for Ak = {A/N) ^ A] = (A/N), the Bethe 
ansatz equations can never yield all eigenstates and eigen- 
values. This becomes clear already on the subspace with 
only one spin fiip. Here the Bethe ansatz equation be- 



A-Nu 



(6) 



which obviously gives only a single solution. 

Therefore integrability ensured by Bethe ansatz breaks 
if all couplings become identical. We now show that in 
this case integrability is ensured by CSCO. In order to 
construct the respective operators we apply the so-called 
binary tree formalism [13 . On the first sight this seems 
to be unnecessary because Gaudin also gave the following 
set of operators which together with the Hamiltonian of 
his central spin model form a CSCO [TO] : 



AiS ' li 



■ E 



AiAj lilj 



Ai 



(7) 



Indeed these operators, which do not play any role con- 
cerning the construction of the eigenstates and eigenval- 
ues of the Gaudin model, obviously become ill-defined in 



Central spins Bath spins 
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FIG. 1: Binary tree with n = A^c + A^ = 2 + 4 = 6 leaves. In 
order to embed the Hamiltonian of an arbitrary central spin 
model with homogeneous couplings into a CSCO using the 
binary tree formalism, two binary trees a, /3 with Nc and N 
leaves respectively must be grafted together. The Hamilto- 
nian results as given in (Isl. 



the homogeneous coupling limit. We restrict ourselves to 
a special case of the binary tree formalism [13 directly 
adapted to our model: Let T be a binary tree with n 
leaves as shown in Fig. [T]for n = 6. A binary tree consists 
of a set of nodes, each of which is connected to exactly two 
following nodes, except for the leaves. If we distinguish 
between a left and a right "child" pi and p2 connected 
to a node p, we arrive at a natural ordering of the leaves. 
We denote the leaves a node p e T is connected to as 
L{p). The node connected to all leaves is called the root, 
denoted by r in the following. Now we associate every 

leaf i with a spin Si and define Hp = (^Xl^^^,^^) Si^ and 



GL(r) 



6'f . It is simple to see that for all p G T 



these operators commute. As every binary tree with n 
leaves has (n — 1) nodes apart from the leaves, we thus 
arrive at exactly n non-trivial, mutually commuting op- 
erators, which indeed form a CSCO. What makes these 
operators "suitable" in the sense of the introduction is the 
fact that they are complete for all spin lengths. In fact 
for any system it is possible to find a CSCO by e.g. con- 
sidering the eigenbasis of the respective Hamiltonian and 
choosing a sufficient number of diagonal matrices with 
only one entry different from zero. We investigated such 
systems for the simple model of two Heisenberg coupled 
spins and found that they consist of more than two oper- 
ators and lose the property of being complete, when the 
spin length is changed. We suppose that sets of commut- 
ing operators can only be complete for any spin length if 
the number of operators is equal to the number of spins. 
Surprisingly, up to our knowledge such a statement has 
not been made so far. 

Now we show how to embed the Hamiltonian of an 
arbitrary central spin model with homogeneous couplings 
in a CSCO. To this end we consider two binary trees 
a and /3 with Nc and TV leaves respectively. Grafting 
them together as shown in Fig. [l] we arrive at a new 
binary tree with Nc^N leaves. If we denote li as 
the Hamiltonian of the associated homogeneous coupling 
model can be written in terms of elements of the CSCO 
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resulting from the binary tree formalism as 

H=^{Hr-Hp)+(^^-^^H^. (8) 

Note that the number of central and bath spins as well 
as their lengths are arbitrary and that there is no further 
restriction to a and P so that indeed there are numerous 
CSCO in which H can be embedded. Furthermore, it 
should be mentioned that by adding JbH^ to ([8| we can 
easily include a homogeneous interaction of strength J5 
between the bath spins. It is simple to find the common 
eigenstates of the respective CSCO [T4]: 

— 5^ I Y\ ' *^P2 5 5 *^P2 I 

i^r,...,^^+2), (9) 

Here Tl = T \ L{r), Tj. = T \ r and Sp denotes the 
quantum number associated with Hp. The complexity for 
calculating the Clebsch-Gordan coefficients is polynomial 
^15 and hence the approach indeed yields integrability. 
The eigenvalues read: 

E{{S,eT,},S:) - ^{Sr{Sr + l)-S^{S^ + l)) 

+ (^-^)^4^. + l) (10) 

Now we relate our above findings to the phenomenon 
of decoherence. The product of two spin operators con- 
sists of flip-flop terms involving ladder operators and a 
coupling of the z-components [16]. In the following we 
evaluate the dynamics for an initial state which is a sim- 
ple product state. In this case all dynamics and hence 
all decoherence is purely due to the flip-flop terms. It 
is well-known that applying a magnetic field B to the 
central spin system strongly suppresses the influence of 
fiip-fiop terms between the central spin system and the 
bath fTP . Here it is usually expected that whenever the 
magnetic field exceeds all other energy scales B ^ l^l, 
a complete neglect of their influence is justified. In the 
following we show that the effect of strong suppression of 
those flip-flop terms actually relies on the inhomogeneity 
of the couplings and is weakened stronger and stronger 
the more couplings are chosen to be equal to each other. 

To this end in Fig. |2] we consider the special case 
Nc = 2 with Si = li = 1/2 and plot the spin dynamics 
for two integrable models (Aj = =: Aj, as explained 
above) with inhomogeneous and homogeneous coupling 
costants. In the first case the coupling constants Aj are 
chosen with respect to a non-uniform distribution so that 
Ai 7^ Aj. For our initial state this case can only be ac- 
cessed via exact diagonalization, strongly restricting the 
size of the system [8 . We therefore illustrate the two sit- 
uations considering a comparatively small system with 
= 2N]j + 1 and Njj = 5. This corresponds to a very 
low bath polarization of 1/N. The initial state of the 
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FIG. 2: Spin dynamics for Nc = 2, N = 2Nd + 1 = 11, where 
Si = li = (1/2), and B = 3.441 A, J = 0.023A. The initial 
state of the system is We consider two 

integrable models fulfilling A\ — A^ —\ Ai with either Ai 7^ 
Aj (upper panel), chosen due to a non- uniform distribution, 
or Ai — Aj (lower panel). Although in both cases B is larger 
than any other energy scale, for homogeneous couplings the 
dynamics still decays, indicating decoherence. 



central spin system is |^-^). We checked the dynamics 
for much larger systems in the homogeneous case using a 
semi-analytical approach based on [16| and did not find 
any qualitative differences. Moreover, non-integrable sys- 
tems with fully inhomogeneous couplings A\. ^ A^, show 
a qualitatively very similar behavior to the integrable 
case of inhomogeneous couplings. A] = =: Ai and 
Ai ^ Aj. Note that all results derived for the special 
case of Nc = 2 and Si = li = 1/2 in the following can 
be directly adapted to the general case of an arbitrary 
number of central spins and arbitrary spin lengths. 

Although the magnetic field is in both cases larger than 
any other energy scale, the dynamics for the inhomoge- 
neous case is completely coherent, whereas in the other 
case it still decays. This means that in the inhomoge- 
neous case the flip-flop terms between the central spin 
system and the bath do not contribute to the dynamics 
in any determinable way. The oscillations are completely 
due to the flip-flop terms between the two central spins. 
A qualitative explanation of the above effect goes as fol- 
lows: Flipping a spin in a magnetic field changes the 
energy E by AE (x B. In order to ensure energy conser- 
vation this change must be compensated. As indicated 
in the upper sketch of Fig. [3) for inhomogeneous cou- 
plings this has to be done by the energy change due to 
the flip of the respective bath spin and the one resulting 
from the central spin flip via the central spin coupling 
term. Hence if the magnetic field exceeds any other en- 
ergy scale, this is impossible and flip-flop processes are 
forbidden by energy conservation (at least in ffist-order 
time-dependent perturbation theory). If we instead con- 
sider homogeneous couplings, this restriction can be cir- 
cumvented by simultaneous flip-flop processes on both of 
the central spins. Here the energy changes due to the 
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FIG. 3: Sketch of flip-flop processes in a A^c = 2 central spin 
system with inhomogeneous (upper panel) and homogeneous 
(lower panel) couplings. For homogeneous couplings, the en- 
ergy changes due to simultaneous central and nuclear spin 
flips can compensate each each other. This is not possible for 
inhomogeneous couplings because the energies of flips diflFer 
for different nuclear spins. 



Of course our results have to be regarded as a first 
indication into this direction and it would be desirable 
to check them for more general external quantities on a 
wider class of systems. As explained above, in ([8| we can 
easily add a term describing an interaction between the 
different bath spins. Hence in an immediate next step 
it would be interesting to check for which types of bath 
terms the Bethe ansatz integrability still holds and if we 
can find effects similar to those described in this paper. 
In this context see e.g. Ref. [6 . 

We thank F. Gohmann for valuable discussions. This 
work was supported by DFG program SFB631. 



central spin flips in the magnetic field and the bath spin 
flips compensate each other as depicted in the second 
sketch of Fig. |3] This is impossible for inhomogeneous 
couplings because the energy change depends on which 
bath spin is flipped. 

This simple effect vanishes for initial states with a fully 
polarized central or nuclear spin system. However, from 
the above explanation it is clear it will still occur if the 
couplings are varied away slightly from complete homo- 
geneity. This means that the more the couplings ap- 
proach the CSCO integrable limit, the less flip-flop terms 
are suppressed by a magnetic field applied to the central 
spin system. This leads to two different interpretations 
of the results, both of which indicate that it is not nec- 
essarily the integrability or non-integrability itself which 
is related to decoherence, as assumed in previous stud- 
ies [6l|7]: (a) As demonstrated above, the influence of a 
magnetic field applied to the central spin system on the 
decoherence properties strongly differs for models which 
are clearly non-integrable or integrable by Bethe ansatz 
and those which are near to the CSCO integrable limit. 
In the first case the dynamics becomes highly coherent, 
whereas in the second case it still decays. This suggests 
that it is the mathematical structure ensuring integra- 
bility, which determines the reaction of a system on an 
external quantity applied to the central system with re- 
spect to the decoherence properties rather than the in- 
tegrability or non- integrability itself, (b) An even more 
general interpretation results from the observation that 
if we apply a magnetic field to the central spin system, 
the non-integrable models as well as those integrable by 
Bethe ansatz keep the respective property, whereas it is 
lost in the CSCO case. Hence the result suggests that if 
a model is close to a limit in which the integrability is 
broken by some external quantity applied to the central 
system, its decoherence properties will be stronger af- 
fected than those of a system near to a limit with stable 
integrability. It is therefore the breaking of integrability 
which has a negative effect on the decoherence properties 
and not the actual integrability or non-integrability. 



[1] V. I. Arnol'd, Mathematical Methods of Classical Me- 
chanics, (Springer, Berlin 1978). 

[2] see e.g. L. D. Faddeev, arXiv:9605187 (1995). 

[3] see e.g. H. Frahm, J. Phys. A: Math. Theor. 26, L473 
(1993). 

[4] R. Steinigeweg, J. Gemmer, and M. Michel, Europhys. 

Lett. 75(3), 406 (2006). 
[5] C. Emary and T. Brandes, Phys. Rev. Lett. 90, 044101 

(2003) 

[6] J. Lages et aL, Phys. Rev. E 72, 026225 (2005); A. Re- 
lano, J. Dukelsky, and R. A. Mohna, Phys. Rev. E 76, 
046223 (2007). 

[7] S. R. Dahmen et al, J. Stat. Mech., P10019 (2004); K. 
Kudo and T. Deguchi, Phys. Rev. B, 68, 052510 (2003); 
R. M. Angelo et a/., Phys. Rev. E 60, 5 (1999); T. Gorin 
and T. H. Seligman, Phys. Lett. A 309, 61-67 (2003). 
[8] J. Schliemann, A. Khaetskii, and D. Loss , J. Phys.: Con- 
dens. Mat. 15, R1809 (2003); W. Zhang et a/., J. Phys.: 
Condens. Mat. 19, 083202 (2007); W. A. Coish and J. 
Baugh, phys. stat. sol. B 246, 2203 (2009); J. M. Taylor 
et al, Phys. Rev. B 76, 035315 (2007). 
[9] F. H. L. Essler et ah, The one- dimensional Hubbard 
model, (Cambridge University Press, 2005). 

[10] M. Gaudin, J. Phys. (Paris) 73, 1087 (1976). 

[11] J. Schliemann, Phys. Rev. B 81, 081301(R) (2010). 

[12] D. Garajeu and A. Kiss, J. Math. Phys. 42, 3497 (2001). 

[13] B. Erbe and H.-J. Schmidt, J. Phys. A: Math. Theor. 43, 
085215 (2010). 

[14] R. Steinigeweg and H.-J. Schmidt, Mathematical 
Physics, Analysis and Geometry 12(1), 19 (2009). 

[15] J. A. De Loera and T. B. McAllister, Exp. Math. 15, 1 
(2006). 

[16] B. Erbe and J. Schliemann, Phys. Rev. B 81, 235324 
(2010). 

[17] see e.g. W. A. Coish, J. Fischer, and D. Loss, Phys. Rev. 
B 81, 165315 (2010). 



